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Abstract 

& . 

Following the construction described in we use the rational map ansatz to construct 

analytically some topologically non-trivial solutions of the generalised SU(3) Skyrme model 

defined by adding a sixth order term to the usual Lagrangian. These solutions are radially 

symmetric and some of them can be interpreted as bound states of Skyrmions. The same 
O ■ 

, ansatz is used to construct low-energy configuration of the SU (N) Skyrme model. 

o 
o 



1 Introduction 

The Skyrme model [0 is widely accepted as an effective theory to describe the low-energy 
properties of nucleons. It was indeed shown || [|, |5j that in the large N c limit, the Skyrme 
model is the low-energy limit of QCD. The classical static solutions of the model describe 
the bound states of nucleons and every configuration is characterised by a topological 
charge which following Skyrme's idea, is interpreted as the baryon charge. 

The Skyrme model can be used to predict the properties of the nucleons within 10 
to 20% M, Q|. To improve these phenomenological predictions various extensions of the 
models have been proposed mostly by adding higher order terms || 0, || || or extra fields 
pl to the Lagrangian. 



The study of the classical solutions of the Skyrme model has been done mostly using 
numerical methods, but recently Houghton et al. [|lTJ showed that the classical solutions of 
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the 577(2) model can be well approximated by using an ansatz that involves the harmonic 
maps from S 2 to S 2 . The harmonic map describes the angular distribution of the solution 
while a profile function describes its radial distribution. This construction was later 
generalised [12| for the SU(N) model using harmonic maps from S 2 to CP 1 ^^ 1 . Moreover, 
it was shown that using a further generalisation of this ansatz one can construct exact 
spherically symmetric solutions of the SU(N) Skyrme model. 

The same method was also used in |13| to construct solutions of another SU(N) 4th 
order Skyrme model. In this paper, we use the same generalised ansatz to construct 
solutions of the sixth order SU(3) Skyrme model and low-energy configurations of the 
SU(N) models defined in fTJ. 

2 The sixth order Skyrme model 

The Skyrme model is described by an SU(N) valued field U(x, t) which, to ensure finite- 
ness of the energy, is required to satisfy the boundary condition U — > I as |x| — > oo, 
where I is the unit matrix. This boundary condition implies that the three dimensional 
Euclidean space on which the model is defined can be compactified into S 3 and as a result, 
the Skyrme field U corresponds to mappings from S 3 into SU(N). As 7T3(SU(N)) = Z 
each configuration is characterised by its winding number, or topological charge, which 
can be obtained explicitly by evaluating the expression 

247T 2 JR3 

where = (<9 M f/)?7 _1 is the right chiral current. Skyrme's ideas was to interpret the 
winding number associated with these topologically non-trivial mappings as the baryon 
charge. 

The generalised sixth order Skyrme model is defined by the Lagrangian 

E = -^-l J dx 3 (^-TrR 2 + l ^Tr[R l ,R J } 2 + ^\Tr[R l ,R (2) 

where this parametrisation of the model is chosen such that A G [0, 1] is a mixing parameter 
between the Skyrme term and the sixth order term: when A = the model reduces to 
the usual pure Skyrme model while for A = 1 the Skyrme term vanishes and the model 
reduces to what we refer to in what follows as the pure Sk6 model. 
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The Euler-Lagrange equations derived from (0) for the static solutions are given by 



di ( Ri — — (1 — A) Rj, [Rj, Ri 



16' 



-A 



Rj, [Rj, Rk][Rk, Ri 



0. 



and the following inequality holds for every configuration 



E > VI - AS. 



(3) 



(4) 



The multi-Skyrmion solutions of the SU(2) Skyrme model have been studied in [14 
where it was shown that they have the same symmetry as the pure Skyrme model. It was 
also shown that the harmonic map ansatz gives a good approximation to the solutions. 



In the next section we describe the harmonic map ansatz. In the third section we 
prove that due to a constraint coming from the sixth order term, the multi-projector 
harmonic map ansatz provides exact solutions only for the SU(3) generalised model. We 
then show that one can nevertheless use the ansatz to construct low-energy configurations 
of the SU(N) models. In the fourth section we look at these configurations for the SU(4:) 
model, while in the last section we look at some special ansatz for the SU(N) model. 



3 Harmonic map ansatz 

The rational map ansatz, introduced by Houghton et al. JTT[] is a generalisation of the 
hedgehog ansatz found by Skyrme 0, to approximate multi-Skyrmion solution of the 
SU(2) model. The ansatz was later generalised by Ioannidou et al. (T| to approximate 
solutions of the SU(N) Skyrme model using harmonic maps from S 2 into CP N ~ l . This 
generalised ansatz is given by 

U(r,d,<p) = e"fWP(o«>)-W 

= e ~ 2 ^ N (I + (e 2 ^ - 1)P(6, ip)) (5) 

where r,6 and (p are the usual polar coordinates. The profile function f(r) must satisfy 
the boundary conditions /(0) = n and lim^oo f{r) = and P{6, ip) is a projector in C N 
which must be a harmonic map from S 2 into CP N ~ l or equivalently a classical solution 
of the 2 dimensional CP 1 ^^ 1 o model. These solutions are well known [|15|, [T(J and to 
construct them it is convenient to introduce the complex coordinate £ = tan(^/2)e* I/: ' 
which corresponds to the stereographic projection of the unit sphere onto the complex 
plane. 
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In these coordinates, P must satisfy the equation 

dP 

and the solutions of that equation are given by any projector of the form 
where h G C N is holomorphic 

dh , . 

M = °' (8) 

The topological charge for the ansatz (|5|), with the prescribed boundary conditions for 
f(r), is given by the winding number of the S 2 — > CP 1 ^^ 1 . This winding number which 
is itself given by the degree of the harmonic function h [13, IB which must then be a 
rational function of £. 

To approximate a solution, one plugs the ansatz into the energy (|2]) and notices 
that if P satisfies (|6]), the integration over the polar angles and the radius decouple. One 
then has to minimise the integral over the polar angles of an expression which depends 
only on P. Taking for P the most general harmonic map of the desired degree, one then 
has to find the parameters of the general map which minimise that integral. Having done 
this, the profile function / is obtained by solving the Euler Lagrange equation derived 
from the effective energy. 

A special case of this construction is the so-called hedgehog ansatz for the SU(2) model 
corresponding to one Skyrmion. In this case, we have h = (1,0* an d after inserting (|7|) 
into (0) the energy reduces to 

E = ^fdr (/> 2 + 2sinV(l + (l-A)/ r 2 ) + (l-A)^ + A^/ r 2 ) (9) 
Sty J r z r z 

and the equation for / is given by 



, , ^ sin'/ , sin 4 A 2 / , sin 4 / 

frr ( y l+2(l-A)^/+A^/J+-/ r (l- A^/ 



+ sin^ / _ _ x + sinV (A + \ = Q _ (1Q) 



This actually corresponds to an exact solution of the model and it is radially symmetric. 
In Figure 1 we present the A dependence of the energy and in Figure 2 we show the profile 
function / and the energy density for the pure Skyrme model, A = 0, and the pure Sk6 



model, A = 1. 
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Figure 1: Total energy of the 1 Skyrmion solution. 




Figure 2: Function profile / and energy density for the 1 Skyrmion solution of the pure Skyrme 
model, A = 0, and the pure Sk6 model, A = 1. . 
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4 Spherically symmetric solutions for the SU(N) model 



In this section we will follow the construction described in 0, to attempt to construct 
solutions of the extended SU(N) Skyrme model using a generalisation of the harmonic 
map ansatz (||). 

To build the new ansatz we need to introduce an operator P + which acts on any 
complex vector u G C N and is defined as 



P + u = d^u — u 
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\u\ 



Taking a holomorphic vector h(£) we then define P^h = h and by induction V& = P+h = 
P + (P^~ 1 h). These N vectors are mutually orthogonal JIB and so the corresponding 
projectors 



P k = P(P k + h) -- 
satisfy the orthogonality relations 



Pjh(Pjhy 
\Plh\ 2 



PkPj 

N-l 
k=0 



- k = 0, 

SijPk 
1 



N-l, 



(12) 



(13) 



as well as other properties discussed in detail in 



The generalised harmonic map ansatz is then defined as 

U = exp{ig (Po-^)+ig 1 (P 1 -^)-...+ig N - 2 (PN-2-jj)} 



-igo/N 



(I + AoPo) e~^ N (I + A 1 P 1 



-igN-2/N 



(I + A N _ 2 P N _ 2 ) (14) 



where gk(r) are AT — 1 profile functions and A k = e %9k — 1. Moreover for the ansatz to be 
well defined, the profile functions ^(r) must be a a multiple of 2ir at the origin and at 
infinity 



To proceed with our construction, it is convenient to rewrite the Euler-Lagrange equa- 
tions of the model (§) using the usual spherical coordinates 



d r lr 2 R r + ^— ^ (A ere + -4^.4 



4 



sin 2 ^^; + 16 A 



sin' 9 



+ B, 



(pip6r8 ) 



sin 9 



d e \ sin 9 



Re H - — [ A r Q r 4 



r 2 sin 9 



:A 



+ 



A 



16 r 2 sin 9 



(B 



rripOtp 



+ 



Btnipr 



ipiprdr , 



+ . o 7, d<pl R v H ; — { A rwr H — -A 6 



sin 2 9 



+ 



A 



16 r- 



{B„.Q V Q + BQQ r 



tpr ) 



(15) 



6 



where Ajy 
that 



Rj, [Ri, Rj] 



and Bjjkik 



Rj, [Rj, Rk] [Ri, Rk] ■ It is fairly easy to show 



N-2 . j 
3=0 V ly 



(16) 



where c/j is the derivative of gj(r) with respect to r. Using the complex coordinates £ and 
£ introduced before we have 



JV-l 



R i = E [e^ - *" - 1 



i=l 



l^-ll 2 



(17) 



and the derivatives with respect to 9 and ip are given by 

1 + KI 2 



ft 



Substituting the above into equations (p~5|) we get 



(9 r 



r 2 P r + (1 - A) 



(I + ICI 



2\2 



(1 + KI 2 ) 2 



+ (P 5 )f + 



(1-A) 



(i + iei 



2\3 



+(1-A) 



(1 + KI 2 ) 2 



HAia-Uui ) + (i-A) 



/ 1 , 



(i + iei 



2\4 



16r 2 



.4 



+ [A-€r]f ) + 



16 



a + iei 2 ) 



2\4 



+ 



ii±Jp!(^[ ( i + iei 2 ) 2 ^- 
(i + iei 2 ) 2 



ft 



2 |^|2 r 2 
(1 + KI 2 ) 2 



4iei 2 



2\2 



+ 



(i+iei 2 ) 2 , 



(1 + l£l ) (B^^ r ^r + B££ r £ r B^ r ^ 



+ 



~^^B^ r ^ 



+ 



ft 



(1 + 161 ) ( — Bg£ r £ r + B^ r ^ r + B^ r ^, 



0. 



In [|T| it is shown that if one takes the special holomorphic vector 

V'o = h = (h , h,..., h N -iY 

where 



(19) 



(20) 



(21) 



and where C^~ l denotes the binomial coefficients, then the terms in ( |T9"[ ) coming from 
the usual Skyrme model, i.e. all the terms except the ones proportional to A/16, are all 
proportional to Pi — Pj_i and p — -L. Using (|13|) one can get rid of the projector Pjv-i and 
(H) will then be the sum of the N — 1 terms Pj — jt for i — ... N — 2, with coefficients 
that depend only on r. This implies that the equations for the Skyrme model reduce to 
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N — 1 ordinary differential equations for the profile functions <?j and their solutions, if 
they exist, will provide us with exact solutions of the SU(N) Skyrme model. 



In what follows we will show that the angular dependence of the terms proportional 
to A in (|T9D, i.e. the terms coming from the sixth order term, is also coming exclusively 
from the projectors Pj — 4? or Pj — Pj_i but that we have to impose an extra constraint 
on the profile functions 

We start by noting that 



— 


to 






* '1 1 


\Vi-i 1 



1=1 

N-l y V \ N-l yyt 

i= i l^-il 2 U. l^-ir 

iV-l y_ yt iV-1 y_ yt 

[i2r, Pf] = 1 E - 9i-lOi) itT 1 12 = E ^ ItV 1 |2 ( 24 ) 

i=l \ v i-l\ i=i 

where Oj = e l ^ 9i ~ 9i - 1 ' — 1. It is then straightforward to check that 

/ IT/. .12 | T/ J2 |y|4 _ |y +i |2 |y|2 \ 



i=l 







\Vi-i\ 2 








\Vi-2\ 2 |1 


^ll 2 



%fr, - %^ = E ( ^ | T r 12It/'|2 + ^j^ji + di ' y 2 ' '/' 2 J (fl - Pi-0 (25) 



where 6j, q and <ij are functions of only. However, as shown in KJ , if Vq is given by 
and (HD then oc (1 + |£| 2 )" 2 and thus 



4 

Furthermore, we have 



3- - % e - rf ) oc (P 4 - P_0 . (26) 



S "^ = i g( e ^ + S ^J^7F (2T) 

with ej = e(^j) and Sj = s(^). But in equation flTP| ) this term appears as 

% [(1 + |£| 2 ) 2 P rr ^-] = 2£ (1 + |e| 2 )S pr€£f + (1 + |£| 2 ) 2 % (B Pr « f ). (28) 

Since % La ^ — 2£ (1 + |£| 2 )~ 3 the only parts of fl28"|) that are non zero are the ones 

V' v ^ - 
that involve the derivatives of ^ '~p with respect to £. Since it can be shown that the 

latter are proportional to E^ 1 C«(l + |£| 2 )~ 2 (P - P_i) where C< = then one sees 

that the term that involves B rr ^^ in (|i~9|) is proportional to (Pj — Pj_i). 

Using similar arguments, it is easy to check that the terms involving B rr ^^, Be^ r t r , 
B^ r ^ r , B^ r ^ r , B^ r ^ r , B^ r ^ r and B^ r ^ r factorise in the same way. 
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There are a few terms in (|T^) which we still have to consider. They involve the 
expressions 



N-l 



i=3 
N-l 



l^ 3 | 2 



(29) 
(30) 



i-3 



where Ki = i (^a, — g>j_iaj). It is clear that these terms will always give a £, £ dependence 
besides the projectors Pi and hence, if we want (|19"|) to reduce to iV — 1 equations that 
involve only the profile functions gi then we have to make sure that (^9[) and fl30|) vanish 
ie we must impose the conditions 



OiKi-iKi-2 — ai_ 2 KiKi-\ 







9i = 9i-2- 



(31) 



This last constraint which is a result of the addition of the sixth order term, implies that 
we can only consider two profile functions go an d g\ and that we should thus have only two 
equations. Unfortunately we have iV — 1 equations which are not compatible with each 
other. From this we see that the ansatz ([5]) will provide exact solutions of the generalised 
Skyrme model for the SU(2) and the SU(3) model only. For larger values of N, the ansatz 
will nevertheless give some low-energy radially symmetric configurations. The SU (2) case 
is nothing but the usual hedgehog ansatz and we will focus on the solutions of the SU(3) 
model in the next section. 



In order to derive the equations for the profile functions, it is convenient to write the 
energy density of the model in terms of (£,£): 

1 



E 



12tt 2 



-(1-A) 



r dr d£d£ Tr 



;i + iei 



2^2 r 



(i + iei 

16r 4 



2\2 



-[R^R^ + X 



64r 4 



2\2 



[R r , , [R r , R^] 



Defining 



as well as Wi = ly.^p (1 — cos(F)) and Wn-i 
above expression can be rewritten as 



|Vjv- 



\V N . 



(N-2 



N-2 



TrR 2 r = ^[T,9i) - E<? 

iV \i=o / i=0 
N-l 

Tr\R^ = -2 £ W { , 

8=1 



2 

i ) 



.(32) 



^ = ft - ft+i for z = 0, . . . , iV - 3, 



Fn-2 — 9N-2 

Wi\ 2 



(33) 



^p- (1 — cos(g)) the terms in the 

(34) 
(35) 
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N-l 



Tr [R rj R&Rr, = -2 £ W* F 2 _ 1; 

k=l 

Tr [Jfc R ( } 2 = 4 (w 2 + " W i+ if + W N _,) , 



(36) 
(37) 



Tr 



7V-2 



4(F 2 iy 2 + J2 {Fi-iWi - FiWi+i) + F 2 _ 2 W 2 N _ X 

(3 



i=l 



In fill it was shown that 



\V k \ 2 
\Vk-i\' 



k(N-k)(l + \Z\ 2 ) 



2\-2 



(39) 



and from this we see that all the terms in ( [32]) are proportional to (1 + |£| 2 ) 2 and that 
after integrating out the angular dependence the energy reduces to 



E= — lr 2 dr< 
6n 



/N-2 



N-2 



N-l 



(1-A) 



N-l 



N 



\i=0 



i=0 
iV-2 



fc=l 



2r 2 



fc=i 



4r 4 
A 



Z 2 + X (Z k 



k=l 



16H 



l) + ^iV-2 f ; 



(40) 



where Z fc = A; (AT — fc)(l — cos(F^_i)). 



In [|J the fields Fi defined by (^3|) were used, and very special solutions were obtained 
by taking F = Fi = . . . = F N _ 2 - It was observed that when Fj(0) = 2tt and Fj(oo) = 
this solution of the SU(N) pure Skyrme model has a topological charge B = ^(iV 2 — 1) 
and has an energy equal exactly to y(A^ 2 — 1) times the energy of the single Skyrmion 
solutions. It is easy to show that, if one uses the same ansatz for the sixth order Skyrme 
model, the profile / = F /2 satisfies the hedgehog profile equation (|i0|) and the energy of 
the configuration is given by E(X) = 4E (X) where E Q (X) is the energy of the hedgehog 
solution for the generalised model. These configurations are not exact solutions, except 
for the 577(3) model. 



To consider the most general ansatz, one can derive from ( [4~0"D the following equations 
for the profile functions Fi, I — 0, ..(N — 2). 



o(] + i) N ~ 2 1 N ~ 2 (i — \) 

- [ ] ^( i + i)F i + 2^^F i + ^^F l (l + l){N-l-l)(l-coBF l ) + 

i=0 k=0 i=k T 



N 
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I -I) 



\ 2 — — \ i — k f r 

~ (I + 1)(N -l-l) sinF, - (I + 1) 2 (N - / - l) 2 (1 - cosF,) sinF, + 

i!_^)(j+l)(jV-J-l) sinF [/(iV-/)(l - cosF,_i) + (/ + 2)(iV-/-2)(l - cosF m )] 

+ A| 2 F(/ + l) 2 (iV - I - 1) 2 (1 - cosF,) 2 - F^Z + l)(iV - l)(N 
(1 - cosF i _ 1 )(l - cosFi) - F l+1 (l + l)(l + 2)(iV - / - 1)(N -1-2) 
(1 - cosF,)(l - cos F m )j + ^ 1 2 F t (l + 1) 2 (N - I - 1) 2 (1 - cosF) 2 - F_ x 
1(1 + l)(N — l)(N — I — 1)(1 - cosF,_i)(l - cosF,) - F m (/ + 1)(Z + 2)(N -l-l) 
(N-l- 2)(1 - cosF,)(l - cosF m )| + ^{2^ + l) 2 (iV - / - 1) 2 (1 - cosF) 
smF t - Pf_ x l(l + 1)(N - l)(N -l-l) sinF^l - cosF) - F 2 +1 (/ + 1)(/ + 2) 
(N-l- 1)(N - I - 2)(1 - cosF) sinF m j = 0. (41) 

When N = 3, the solution of the 2 equations lead to exact solutions of the model, 
while for larger values of N, the ansatz (|L4"D corresponds to low-energy configurations. 

We would like to point out at this stage that as proved in []TJ, the topological charge 
for the configuration ([14]) is given by 

N-2 

B= J2^k(Fi- sin F^ (42) 

i=0 

where 

V* = I (43) 



'4tt 2 J \P£h\ 2 

takes integer values given by the degree in £ of the wedge product [0 of h and its 
derivatives 

V k = ^-deg(h (k) ) h {k) = hAd+hA... Ad k + h k = 0,..N-l. (44) 

2n 

Each configuration is thus characterised by the boundary conditions for the profile 
function Fj and we can without loss of generality impose the condition lirrir^ooF^r) = 0. 
For the configuration to be well-defined at the origin we must also impose a condition of 
the type 

Fi(0) = ni 27T (45) 

where the m G N. 



11 



5 Radially symmetric 577(3) Solutions 

To describe the solution of the 577(3) model, we use the profile F = F and g — F% and 
the energy fllPf) simplifies to 



E = ^Jr 2 drl^(g 2 + F 2 + gF) + ^((l-cosF)((l-\)F 2 + 4) 



2 

+ (1 - cos0)((l - A)# 2 + 4)) + (1 - A) — ((1 - cosF) 2 - (1 - cosF)(l - cose/) 

+ (1 - cosg) 2 ) + A (F 2 (1 - cosF) 2 + g 2 (l - cos^) 2 - (1 - cosF)(l - cosg)gF) |. 

(46) 

The equations for the profile function F and g are then given by 

9rr + \f„ + y + 2 j + ^((1 - A)(l - cos g)g rr + ^ sin#((l - \)g 2 - 4) 
+ l - sin g((l - X)g 2 - 4) + (1 - A) - cosF) - 2(1 - cos g)) sin(^) 



3A / 

■g^(l - cos5f)f 2(sin5(5( 2 + (1 - cosg)(g rr - 2-j- 



- sin FF 2 - (1 - cos F)(F rr - 2— ) ] = (47) 



r 



F rr + - # rr + 2 — + ^ + — sin F((l - A)F - 4) + 2(1 - A) (1 - COS F)F r , 
2 f f 4r ^ 

3 3A / 

-(1 - A) — (2(1 - cosF) - (1 -cos#))sinF + — (1 -cosF) 12(sinFF 2 
2r 4 8r 4 y 

+ (l-cosF)(F rr -2^))-sin^ 2 -(l-coss)(2 rr -2^))) =0. (48) 

r r J 

The topological charge of the solution now reads 

B = i ( (F - sinF) |™ + (<? - sin(^))|^ (49) 
and if we take the boundary conditions 

F(0) = n F 2ix 

9(0) = n g 2n (50) 

where rip and n g are integers, we have B = 2(nj + n g ). When np and n g are of opposite 
signs, we can interpret the solutions as a mixture of Skyrmions and anti Skyrmions. 
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In Table 1, we give the energy of the hedgehog solution (B = 1) for the SU(2) model. 
This solution is an embedded solution of any SU(N) model and it is the solution with 
the lowest energy. We thus use it as the reference energy for all the other solutions. 

In Table 2 we present the properties of the different solutions for the 577(3) models. 
The first two columns specify the boundary condition of the solution, and the third 
columns gives the topological charge of that solution. In column 4 and 5 we give the 
energy of the solutions for the pure Skyrme model and the pure Sk6 model while column 
6 and 7 give the corresponding relative energy per Skyrmion, that is the energy divided by 
the energy of the single Skyrmion and the total number of Skyrmions. For the solutions 
corresponding to the superposition of Skyrmions and anti-Skyrmion, we define the total 
number of Skyrmions as the total number of Skyrmions and anti-Skyrmions. Notice that 
the cases n g = 0, np = 1 and n g — 1, np — correspond to the same solution modulo an 
internal rotation. 

In Figure 3, we present the energy of the 3 different types of solution as a function of 

A. 

5 
4.5 

4 
3.5 

w 

3 
2.5 
2 

0.2 0.4 0.6 0.8 1 

X 

Figure 3: Energy of the SU(3) solution for the boundary conditions (A) np = 0,n g = 1, (B) 
np = l,n g = 0, (C) np = l,n g = —1, (D) np = l,n g = 1. 

6 Low Energy SU(4:) Configurations 

As was shown in the last two sections, the ansatz (TTJ|) provides an exact solution of 
the sixth order model only for the SU(3) model, or when A = 0, that is for the usual 
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SU{2) 


Energy 


n g 


B 


E(0) 


E(l) 


1 


1 


1.2315 


0.9395 



Table 1: Topological charge and Energy of the hedgehog SU(2) solution. 



SU{3) 


Total Energy 


Relative 


Energy 


rip 


n g 


B 


E(0) 


E(l) 


£? fl (0)/(|B|^(0)) 


E B {l)/{\B\Ei{l)) 


1 




1 


4 


4.928 


3.758 


1 


1 


1 







2 


2.377 


1.819 


0.965 


0.968 







1 


2 


2.377 


1.819 


0.965 


0.968 


1 




-1 


2-2 


3.862 


3.191 


0.784 


0.849 



Table 2: Topological charge and Energy of some SU(3) solutions. 

Skyrme model. For the SU(N) model with N > 4, the ansatz still produces low-energy 
configurations. In particular, when A is small, we can expect the ansatz to be very close to 
an exact solution. In this section we look at some configurations of the SU{4) model. For 
this model, we have three profile functions F , Fx and F 2 and the energy for the general 
ansatz (|i~4"D is explicitly given by 



E 



= Jr 2 dr^ (3F 2 + 4F 2 + 3F| + 4F A + 4FxF 2 + 2F F 2 ) + 1 [3(1- cos F ) 
+4(l-cosF 1 ) + 3(l-cosF 2 )] + (1 - A) [3F 2 (l-cosF )+4Fftl-cosFi}f 
3F 2 2 (l-cosF 2 )l + — r{9 (1 - cosF ) 2 + 16 (1 - cosFx) 2 + 9 (1 - cosF 2 ) 2 
-12(1 - cosF )(l - cosFx) - 12(1 - cosPi)(l - cosF 2 )}| 

{ g ^o (! - cosF ) 2 + 16F 2 (1 - cos^) 2 + 9F|(1 - cosF 2 ) 2 
-12^1(1 - cosF )(l - cosFx) - 12FxF 2 (l - cosFiXl - cosF 2 )} 1 (51) 



from which we can derive the following equations 
/3A(l-cosF ) 2 2(l-A)(l-cosF ) 



2r 4 



+ 



1 F n 



'2 A(l - cosF )(l - cosFi) 



+ -F 2 
3 



4 sin F 6 F + AF\ + 2F 2 (1_- A)F 2 sin F 

nn2 7" 
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+ (1 - A)^(4(l - cosF) - 6(1 - cosF )) + A (1 7^ ( 3 -F*sinF - F^smF, 
-A (1 ~^ sFo) (3F (1 - cosFq) -2F 1 (1 - cosF 1 )) = 0, (5: 



(I 3A(l-cosF )(l-cos(F 1 ) \ tf , / 2A(l-cosF 1 ) 2 2(1 - A)(l - cosfi) U 

/l 3A(1 -cosFi)(l -cosF 2 )\ p (1 - A)F 2 sinFi F + 2F 1 + F 2 4 sin Fx 
\2 4r 4 y 2 r 2 r r : 

+ (1 - A)^i (3(1 - cosFq) + 3(1 - cosF 2 ) - 8(1 - cos Fx)) 
A 3 • „ x 3 



2 



(1 - cos F) (4Fx(l - cos F) - - F (l - cosF ) - - F 2 (l - cosF 2 )) 



+^(1 - cos Fi) ^2 F 2 sin F - ^Fq sin F - ^F 2 2 sin F 2 ) = (53) 

and 



'2 A(l-cosF 1 )(l-cosF 2 )\ /3A(1 - cosF 2 ) 2 2(1 - A)(l - cosF 2 ) \ 

^3 H J Fl+ i 2^ + V- + 1 J 2 

iF + 2F ° + 4F + 6F 2 _ 4 sn^ + (1 - A)F| smF 2 + _ srn^ _ ^ 
3 3r H r 4 



—6(1 — cosF 2 )) — A— — ^31 (3F 2 (1 -cosF 2 ) -2F(1 -cosFi)) 

+ A (1 ~7 F2) (? F 2 sin F 2 - F 2 sin F x ) = 0. (54) 

Describing the boundary condition for the profile functions as before, Fj(0) = nj27r, the 
topological charge is given by 

B = 3n + 4m + 3n 2 . (55) 



In Table 3 we present the energy values of various types of configurations when A = 
and A = 1. We notice that when A = 0, the solutions are symmetric under the exchange 
fo /2, but that the sixth order term breaks the symmetry. This results in a difference of 
energy between the configuration with n — 0, rt\ — 0, n2 — 1 and n = 1, ni = 0, n2 = 
as well as between the configurations with n — 1, rt\ — 1, n2 — and n = 0, ni = 1, n2 = 
1. In Figure 4, we present the curve for the energy of the configurations as a function of 
A. 



15 



SU{4) 


Total Energy 


Relative Energy 


n 


ni 


n 2 


B 


E(0) 


E(l) 


^(ova^i^o)) 


E B (i)/(\B\m)) 








1 


3 


3.51739 


2.66653 


0.95210 


0.94598 


1 








3 


3.51739 


2.72915 


0.95210 


0.96819 





1 





4 


4.78807 


6.33322 


0.97204 


1.68507 


1 





1 


6 


7.22464 


6.04604 


0.97780 


1.07244 


1 


1 





7 


8.45219 


6.62998 


0.98052 


1.00802 





1 


1 


7 


8.45219 


7.28058 


0.98052 


1.10694 


1 


1 


1 


10 


12.311 


9.39605 


1 


1 



Table 3: Topological charge and Energy of some SU(4) configurations. 




Figure 4: Energy density of the SU(4) multi-projector ansatz (a) no = 0, n\ = 0,ri2 = 1; 
(b) n = l,ni = 0, n 2 = 0; (c) n = 0, n\ = l,n 2 = 0; (d) n = l,ni = 0, n 2 = 1; (e) 
n = l,ni = l,n 2 = 0; (f) n = 0, n\ = l,n 2 = 1. 
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7 SU(N) Low Energy configuration 



After inserting the ansatz (|5|) in the full equation for the SU(N) model, we found that 
we had only two independent profile functions go and g\ and that the ansatz would only 
provide solutions for the SU(3) model. One can nevertheless use the SU(N) ansatz 
to compute low energy configurations. For example if we consider the reduced ansatz 
defined by (|) together with the constraint gi = g i+2 and define the profiles F = go — g± 
and g = #jv-2 we can minimise the energy (HDI) and solve the equations for F and g for 
various boundary conditions. We found that to get configurations corresponding to a 
bound state, i.e. a configuration with an energy per Skyrmion smaller than the energy of 
the hedgehog solution, we must take rip = and n g = 1. The energies that we found are 
given in Table 3. 





Total Energy 


Relative Energy 


Model 


B 


E(0) 


E(l) 


E B (G)l(\B\Ei{Q)) 


E B {l)/{\B\Ei{D) 


SU(3) 


2 


2.377 


1.819 


0.965 


0.968 


SU(A) 


3 


3.624 


2.759 


0.981 


0.979 


SU(5) 


4 


4.811 


3.632 


0.977 


0.966 


SU(6) 


5 


6.015 


4.518 


0.977 


0.962 



Table 4: Topological charge and energy for the reduced ansatz with np = and n„ = 1. 



In Figures 5 and 6 we present the profile and the energy density for different values 
of N and for A = 0.5. It shows that the energy density has the shape of a hollow sphere 
of radius r = 0.7\/N. The profile g has the same shape for all values of iV but is shifted 
to the right as N increases. The profile F on the other hand is also shifted as the shell 
radius increases, but its amplitude decreases like l/N 2 ; note that in Figure 6, the profile 
for iV = 100 and iV = 200 have been multiplied by 100 to make them visible. For other 
values of A the graphics look very much the same except that the shell radius and width 
are slightly different, but the conclusions remain the same. 

Figure 6.b suggests to simplify the ansatz further for large iV by taking F(r) = 0. 
This implies that g^ = g Vi and the multi-projector ansatz (|5|) becomes 

U = exp{-ig(P N . 1 -I/N)) (56) 
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5 10 
r 



15 




(a) (b) 

Figure 6: Profile (a) g and (b) F of the multi-projector solution with np = 0, n g = 1, A = 0.5. 
(A) F for N=10, (B) F for N=20, (C) F for N=50, (D) 100 x F for N=100, (E) 100 x F for 
N=200. 
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where Pn-i can also be written as 

hh) , . 

p N _ x = 57 

where h is equal, up to a unitary rotation, to the complex conjugate of the holomorphic 
vector V Q defined in ([^H]) : h = AV for some A e SU(N) with = <%A = 0. This is 
shown by using the fact that Pn-i is an anti-holomorphic projector and that solving 
(ft9p recursively we have 

n4i'= ( ^r-iV i+itir "'"" (58) 

and so |Vat_i| 2 = (JV — 1)! 2 |1 + |£| 2 | 1-Ar - Knowing that up to an overall coefficient | Vjv__ 1 1 2 
is a polynomial in £ of degree N — 1, we can conclude that up to a unitary iso-rotation, 
Vjv_i is equal to the complex conjugate of Vo- 

The topological charge of the anti-holomorphic projector Pn-i is equal to 1 — iV and 
as the profile function is —g, the baryon number for this configuration is AT — 1. The 
ansatz Q55"D is not a solution, but its energy 

E = ^/r 2 rfr|^^^ + ^ + (iV-l)(l-cos^)((l-A)^ + 4) 

- WO- - cos ^) 2 ((! " A ) + i^ 2 ) } (59) 
can easily be computed by solving the equation 



./V 1 

2g rr + 4^ + — ((1 - A)(l - cos g)g rr + -ang({l - X)g 2 r - < 
+ ^- 4 N(N - 1)(1 - cos^)(sin^ r 2 + (1 - cosg)(g rr - 2^)) = 0. (60) 



In Figure 7, we present the relative energy, E(X) / (Eb=i(X)(N — 1)), of this configura- 
tion as a function of for different values of A. We see that this configuration corresponds 
to a bound state of Skyrmions and that the energy per Skyrmion decreases with N. The 
energy of this configuration corresponds to an upper bound for the energy of the B = N—l 
radially symmetric solution of the SU (N) model and these configurations correspond to 
bound states of Skyrmions for all values of A^ and all values of A. As every SU (p) solution 
can be trivially embedded in an SU(q) solution when p < q we can claim that for every 
B < N the SU(N) model has a radially symmetric solution of charge B corresponding to 
a bound state. With the exception of the hedgehog solutions, these solutions are expected 
to be unstable when the radial symmetry is broken as their energies are larger than the 



known SU{2) solutions [jT4 
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0.8 

w 

0.7 
0.6 
0.5 

50 100 150 200 

N 

Figure 7: Energy E/(E B=1 (N -1)) of the SU(N), configuration @ for (a) A = 0, (b) A = 0.25, 
(c) A = 0.5, (d) A = 0.75, (e) A = 1 

8 Conclusions 

In this paper we have shown how to construct some radially symmetric solutions of the 
SU(3) sixth order Skyrme model. The construction is similar to the one used for the 
pure Skyrme model in except that, because of an extra constraint, the construction 
only works for the SU(3) model. The same ansatz can nevertheless be used to compute 
low-energy configurations of the SU(N) model. In particular we showed that for every 
N there is a radially symmetric solution of charge B < N which corresponds to a bound 
state of Skyrmion. 
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